Random-matrix theory of Andreev reflection from a topological superconductor 



C. W. J. Beenakker, J. P. Dahlhaus, M. Wimmer, and A. R. Akhmerov 
Instituut-Lorentz, Umversiteit Leiden, P. O. Box 9506, 2300 RA Leiden, The Netherlands 

(Dated: December 2010) 

We calculate the probability distribution of the Andreev reflection eigenvalues _R„ at the Fermi 
level in the circular ensemble of random-matrix theory. Without spin-rotation symmetry, the statis- 
tics of the electrical conductance G depends on the topological quantum number Q of the supercon- 
ductor. We show that this dependence is nonperturbative in the number N of scattering channels, 
by proving that the p-th cumulant of G is independent of Q for p < N/d (with d = 2 or d = 1 in the 
presence or in the absence of time-reversal symmetry). A large- A effect such as weak localization 
cannot, therefore, probe the topological quantum number. For small N we calculate the full dis- 
tribution P{G) of the conductance and find qualitative differences in the topologically trivial and 
nontrivial phases. 
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I. INTRODUCTION 

The random-matrix theory (RMT) of quantum trans- 
port studies the statistical distribution of phase coher- 
ent transport properties in an ensemble of random scat- 
tering matrices. The theory finds a major application 
in the context of chaotic scattering, because then the 
ensemble is fully specified by fundamental symmetries 
— without requiring microscopic input. ^ Since scatter- 
ing phase shifts for chaotic scattering are uniformly dis- 
tributed on the unit circle, such ensembles are called 
"circular", following Dyson who first introduced these 
ensembles in a mathematical context.^ The circular en- 
sembles have been successful in describing experiments 
on low-temperature electrical and thermal conduction in 
quantum dots, which are confined geometries connected 
by point contacts to metallic or superconducting elec- 
trodes. For a recent overview of the field we refer to 
several chapters of a forthcoming Handbook.'^ 

While metallic quantum dots are characterized by the 
three circular ensembles introduced originally by Dyson, ^ 
superconducting quantum dots are described by four dif- 
ferent ensembles discovered by Altland and Zirnbauer<^ 
The classification of the superconducting ensembles is 
based on the presence or absence of time-reversal and 
spin-rotation symmetry, as summarized in Table U] The 
symmetry classes are called D, Dili, C, and CI, in a no- 
tation which originates from differential geometry^i The 
corresponding circular ensembles, in the nomenclature of 
Ref. are the circular real (CRE) and circular quater- 
nion (CQE) ensembles in class D and C, respectively, and 
their time-reversal invariant restrictions (T-CRE and T- 
CQE) in class Dili and CI. 

In a remarkable recent development, it was found 
that some of these symmetry classes divide into dis- 
junct subclasses, characterized by a topological quantum 
numberi^""— For a quantum dot connected to a supercon- 
ducting wire, such as shown in Fig. [U the topological 
quantum number Q = —1 or -1-1 depending on whether 
or not the quantum dot has a bound state at zero ex- 
citation energy. (The state is only quasi-bound if the 




FIG. 1: Quantum dot with an interface between a normal 
metal (N) and a superconductor (S, shaded region). Andreev 
reflection at the NS interface converts a normal current (car- 
ried by electron and hole excitations e and h) into a super- 
current (carried by Cooper pairs). The conductance G is the 
ratio of the current / into the grounded superconductor and 
the voltage V applied to the quantum dot via an A-mode 
point contact to a normal metal electrode (narrow opening 
at the left). The system is in a topologically nontrivial state 
if it supports a quasi-bound state at the Fermi level. This 
is possible if spin-rotation symmetry is broken by spin-orbit 
coupling. In the configuration shown in the figure (with a sin- 
gle NS interface), time-reversal symmetry should be broken 
to prevent the opening of an excitation gap in the quantum 
dot. In the presence of time-reversal symmetry a second NS 
interface, with a tt phase difference, can be used to close the 
gap. 



quantum dot is also connected by a point contact to a 
normal-metal electrode.) Because of particle-hole sym- 
metry, such a fermionic excitation is equal to its own 
antiparticle, so it is a Majorana fermion. There is now 
an active search for the Majorana bound states predicted 
to appear at the ends of superconducting wires without 
spin-rotation symmetryiiSr— 

The RMT of the four superconducting circular en- 
sembles was developed for the quasiparticle transmission 
eigenvalues in Ref. [H, and applied to the thermal con- 
ductance. The probability distribution of this transport 
property does not depend on the topological quantum 
number, so it was not needed in that study to distin- 
guish the topologically nontrivial ensemble (with a Majo- 
rana bound state) from the topologically trivial ensemble 
(without such a state). 
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Symmetry class 


D 


Dili 


C 


CI 


Ensemble name 


ORE 


T-CRE 


CQE 


T-CQE 


Particle-hole symmetry 






Time-reversal symmetry 


X 


T Sy 


X 


T 

r — r 


Spin-rotation symmetry 


X 


X 


/ 


/ 


topological quantum number Q 


Detr 


Pf iEyr 


X 


X 


degeneracy of iin 7^ 0, 1 


2 


2 


2 


2 


degeneracy of -Rn = 0, 1 


1 


2 


2 


2 



TABLE I: Classification of the symmetries of the reflection matrix r for a normal-metal-superconductor junction. See Sec. |II] 
for explanations. 



The electrical conductance G, in contrast, can probe 
the presence or absence of the Majorana bound state 
through resonant Andreev reflectionJ^iil This applies 
also to a quasi-bound state,— so even if the A^-mode 
point contact to the normal metal has a conductance 
which is not small compared to /h — although the 
effect of Q on G diminishes with increasing N . As we 
will show in this paper, the Q-dependence of the con- 
ductance distribution P(G) in the circular ensembles is 
nonperturbative in N: Cumulants of order p are identi- 
cal in the topologically trivial and nontrivial phases for 
p < N {N/2) in the absence (presence) of time-reversal 
symmetry. 

We derive this nonperturbative result by exactly com- 
puting (in Sec. Illll and App. IXj) the joint probability dis- 
tribution P{{Rn}) of the Andreev reflection eigenvalues 
i?2, • ■ • Rn (eigenvalues of the product rj^^r^^^ of the 
matrix r^e of Andreev reflection amplitudes). This prob- 
ability distribution was only known previously for one 
of the ensembles (CQE) without a topological phase^^P 
while here we calculate it for all four superconducting en- 
sembles, including the possibility of a topologically non- 
trivial phase. 

The distribution of the electrical conductance G oc 
J2n follows upon integration over the i?„'s (Sec. IIVI) . 
For small N we obtain a closed-form expression for P{G) 
in the two cases Q — ±1, and we demonstrate that, upon 
increasing A^, first the average becomes the same, fol- 
lowed by variance, skewness, kurtosis, etc. A proof for 
arbitrary A^ of the Q- independence of low-order cumu- 
lants is given in App. [BJ 

We conclude in Sec. |V] with a comparison of the analyt- 
ical predictions from RMT with a computer simulation of 
a system that is of current experimental interest (a disor- 
dered semiconductor wire on an s-wave superconducting 
substrate, with Rashba spin-orbit coupling and Zeeman 
spin splitting) ji^iii 



II. ANDREEV REFLECTION EIGENVALUES 

Andreev reflection of electrons injected through a point 
contact towards a normal-metal-superconductor (NS) in- 
terface is described by the N x N subblock r^e of the 



2N X 2N reflection matrix r, 

\rhe rhhj 

The labels e and h refer to electron and hole excitations 
in the normal-metal point contact, each of which can be 
in one of A^ modes. We allow for spin-orbit coupling, so 
N counts both orbital and spin degrees of freedom. The 
Andreev reflection eigenvalues i?„ {n — 1,2,...N) are 
the eigenvalues of the Hermitian matrix product r^^^r^^^. 

At excitation energies e below the superconducting gap 
A there is no transmission of excitations through the su- 
perconductor (assuming that it is large compared to the 
penetration depth ^0 — hvp/^)- The reflection matrix 
is then a unitary matrix, — r^^ . An electrical cur- 
rent / can still enter into the superconductor, driven by 
a voltage difference V with the normal metal. The elec- 
trical conductance G = I/V is fully determined by the 
Andreev reflection eigenvalues,— 

N 

G/G^ =N-Tr rlr,, + Tr rlj^^ = 2 ^ i?„. (2.2) 

n=l 

The conductance quantum is Go = /h and the factor 
of two accounts for the fact that charge is added to the 
superconductor as Cooper pairs of charge 2e. (Spin is 
counted in the sum over n.) 

The relation (|2.2I) holds in the limit of zero voltage and 
zero temperature, when the reflection matrix is evaluated 
at the Fermi level (e — )■ 0). The subblocks in Eq. (|2.ip 
are then related by particle-hole symmetry. 

Time-reversal symmetry imposes a further constraint on 
the reflection matrix, 

r = S^r^S^, (2.4) 

with Yiy = ay(BiJy(B- ■ - (BcTy and Cj, a 2 x 2 Pauli matrix. 
(The superscript T denotes the transpose.) 

The Andreev reflection eigenvalues are all twofold de- 
generate in the presence of time-reversal symmetry. This 
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is the usual Kramers degeneracy, which follows directly 
from the fact that T,yr is an antisymmetric matrix [T,yr — 
-(Eyr)"^] when Eq. ^ holds.-^ 

Remarkably enough, a twofold degeneracy persists 
even if time-reversal symmetry is broken. More precisely, 
as discovered by Beri,— if i?„ is not degenerate then it 
must equal or 1 . This follows from the antisymmetry of 
the matrix ?'J[^''ee, which is required by particle-hole sym- 
metry and unitarity irrespective of whether time-reversal 
symmetry is present or not. The degeneracy of the An- 
dreev reflection eigenvalues i?„ ^ 0,1 is remarkable be- 
cause the eigenvalues of the Hamiltonian are not degen- 
erate in the absence of time-reversal symmetry. To dis- 
tinguish it from the Kramers degeneracy, we propose the 
name "Beri degeneracy" . 

The determinant of r is real due to particle-hole sym- 
metry, and therefore equal to either +1 or —1 due to 
unitarity. The topological quantum number Q = Dot r 
distinguishes the topologically trivial phase {Q — 1) from 
the topologically nontrivial phase {Q = —1)^^^ This 
quantum number is ineffective in the presence of time- 
reversal symmetry, when Kramers degeneracy enforces 
Detr = 1. In that case the Pfafhan (the square root of 
the determinant of an antisymmetric matrix) can be used 
instead of the determinant to identify the topologically 
nontrivial phase;^ Q — Pi iYiyV equals -1-1 or —1 depend- 
ing on whether the superconductor is topologically trivial 
or not. 

A topologically nontrivial superconductor has a (pos- 
sibly degenerate) bound state at £ = 0, consisting of an 
equal-weight superposition of electrons and holes from 
the same spin band. It is the tt phase shift upon reflection 
from such a Majorana bound state which is responsible 
for the minus sign in the topological quantum number;^- 

These properties of the reflection matrix in the ab- 
sence of spin-rotation symmetry are summarized in Ta- 
ble HI For completeness, we also include in that table the 
case when there is no spin-orbit coupling. In that case 
it is sufficient to consider only the orbital degree of free- 
dom, with a two- fold spin degeneracy of all i?„'s. The 
conditions for particle-hole symmetry and time-reversal 
symmetry are then both different from Eqs. (|2.3p and 
p.4p . given respectively by 



"^hh — Ke^ ''"eh — (2-5) 

r = r^. (2.6) 

As a consequence, the determinant of r is now always 
+1, while the Pfaffian does not exist (for want of an an- 
tisymmetric matrix). Without spin-orbit coupling An- 
dreev reflection exclusively couples electrons and holes 
from opposite spin bands, which prevents the formation 
of a Majorana bound state at the NS interface. 



trical conduction in the next section. For each symme- 
try class we first determine the polar decomposition of 
r in terms of the R^s. The resulting invariant measure 
d^{r) cx P{{Rn}) n,i dRn then gives the probability dis- 
tribution P{{Rn}) of the Andreev reflection eigenvalues 
in the corresponding circular ensemble. 



A. Class D, ensemble CRE 

In the absence of time-reversal and spin-rotation sym- 
metry the scattering matrix has the polar decomposition 



U 

U* 




(3.1) 



The N X N matrices U, V are unitary and the N x N 
matrices A and F are real, to satisfy the particle- hole 
symmetry condition (j2.3p . Unitarity of r requires, in ad- 
dition to F'^F + A^A = 1, that A'^F = -F'^A is antisym- 
metric. As derived in Ref. [2^ the matrices A and F must 
therefore have a 2 x 2 block diagonal structure. 

For N = 2M even and Q = 1 one has A = Am, F = Fm 
with 

M / . \ M 

AM = 0rY" . U0^osina„, (3.2) 

Af / \ M 

Fm = ^ ^°^""U0za,cosa„. (3.3) 
^ \-cosa„ / ^ 

n—l \ / n— 1 

The 2x2 Pauli matrices are cr^, ct^, az (with (Tq the 2x2 
unit matrix). The real angles a„ € (0, 27r) determine the 
Andreev reflection eigenvalues i?„ = sin^ . These are 
all twofold degenerate. 

The parameterization derived in Ref. [1^ has Am oc iay 
and Tm oc cto- The present, equivalent, form is chosen 
because it is more easily extended to symmetry class Dili 
(where an additional symmetry condition applies). For 
the same reason, we parameterize the diagonal entries in 
terms of the angles q;„ rather than in terms of y/Rn and 
Vl — Rn- (The sign of the terms sina„, cosa„ cannot 
be fixed in class Dili.) 

To check that the polar decomposition p.l|) - p.3|) in- 
deed gives Det r = 1 , one can use the identities Det AB = 
(Det^)(DetS) and 



Det 1^ =I)ct{AD- ACA-^B) 



For N — 2M even and Q = —1 one has 



(3.4) 



III. RANDOM-MATRIX THEORY 



A = Aa/_i ® 





,0 1, 



F = Fa/_i ® 



1 0' 
,0 0, 



(3.5) 



In this section we calculate the distribution of the An- 
dreev refiection eigenvalues, which we then apply to elec- 



so in addition to M — 1 two-fold degenerate eigenvalues 
Ri, R2, ■ ■ ■ Rm-1 there is one nondegenerate eigenvalue 
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equal to and one nondegenerate eigenvalue equal to 
1. It is this unit Andreev reflection eigenvalue which 
contributes a factor —1 to Detr. 

For N — 2M + 1 odd there are M two-fold degen- 
erate eigenvalues i?2, ■ • • Rm plus one nondegenerate 
eigenvalue equal to<7=(l — Q)/2, 



A = Am ® (g), r = Tm ® (1 - q). 



(3.6) 



The nondegenerate eigenvalue equals 1 in the topologi- 
cally nontrivial phase and otherwise. Again, it is the 
unit Andreev reflection eigenvalue which gives Det r = 
-1. 

The calculation of the invariant measure from these 
polar decompositions is outlined in App. El The result- 
ing probability distributions of the twofold degenerate 
Andreev reflection eigenvalues in the CRE are 

M 

P{{Rn}) OC \{{R^- Rj)\ 
i<] = l 

if iV = 2M and Q = 1, (3.7) 

p{{Rn})^ n (R-Pi)' n Rii^-Rk)', 

i<j = l k=l 

if iV = 2M and Q = -1, (3.8) 

M M 

P{{Rn})^ n {R^~R,yll^Rk~U^~Q)]^ 

i<j = l k=l 

if iV = 2A/ + 1. (3.9) 

The degenerate Andreev reflection eigenvalues repel each 
other with the fourth power of their separation. In addi- 
tion there is a repulsion with the second power of the sep- 
aration to each of the nondegenerate eigenvalues, pinned 
at or 1. 



not an additional degeneracy. The topological quantum 
number^^ Q — Pi iTiyV can be calculated using the iden- 
tity 



Pf 



/ a b c\ 

—a Ode 
-b -d f 
V-c -/ 0/ 



af — he + cd, 



(3.11) 



for scalars a,b,c,d,e, f, and also the formulas 

pfxrxT = (Detx)(Pfy), Pf©„r„ = n„Pf^« 

(valid for arbitrary square matrices X and antisymmet- 
ric matrices Y,Yn). 

For Q = 1 we take A — Am and F — Tm from 
Eqs. (1321) and (1331). The Pfafhan of the matrix (|XTU1) 
is always -1-1, so this polar decomposition describes the 
topologically trivial phase. In order to reach the topo- 
logically nontrivial phase, we include a twofold degen- 
erate eigenvalue equal to unity, but with a az matrix 
rather than a cfq matrix: A = Am-i ® diag(l,— 1), 
r = Tm-1 ® diag (0, 0). Then the Pfafhan is -1. 

As derived in App. \^ the distribution of the M de- 
generate Andreev reflection eigenvalues in the T-CRE 
is given most compactly in terms of the variables ^„ = 
sina„ e (—1, 1) (with i?„ = f^). For Q = I the result is 



M 



Pil^n}) CX Yl (e. 



(3.12) 



Notice that there is no repulsion of pairs of Andreev re- 
flection eigenvalues: If — > — then Ri — > Rj and yet 
the probability distribution does not vanish. 

For Q = — 1 one pair of eigenvalues is pinned at Rm — 
1 =^ Cm = 1- The distribution of the remaining M — 1 
degenerate eigenvalues is 



M-l 



M-1 



B. Class Dili, ensemble T-CRE 

In the presence of time-reversal symmetry the scatter- 
ing matrix should also satisfy the condition (|2.4|) . which 
implies that iSj,r is antisymmetric. The polar decompo- 
sition which respects this symmetry condition (as well as 
the condition (|2.3p for particle-hole symmetry) is 



iTiyr 



n 

, n 




(3.10) 



with fl SLii N X N unitary matrix. Unitarity still requires 
that A"'"r is antisymmetric, while time-reversal symmetry 
requires F"^ = -F, A""" = A. 

The number of channels N — 2M is even, with AI 
the number of channels per spin. Each reflection eigen- 
value has a two-fold degeneracy, including those equal to 
or 1. This Kramers degeneracy due to time-reversal 
symmetry2^ in class Dili replaces the Beri degeneracy 
due to electron-hole symmetry^'^ in class D — it is 



(3.13) 



i<j=i 



fc=i 



While pairs of Andreev reflection eigenvalues i?„ G (0, 1) 
in the T-CRE do not repel each other, they are repelled 
from the eigenvalue pinned at Rm = 1, with the same 
second power of the separation as in the CRE. 



C. Class C, ensemble CQE 

For completeness we also consider the two symmetry 
classes C and CI without a topological phase. Then spin- 
rotation symmetry is preserved, so it is sufficient to con- 
sider a single spin degree of freedom, say an electron in 
the spin-up band coupled to a hole in the spin-down band. 
The reflection matrix for this scattering process has di- 
mension 2M X 2M, where M only counts the orbital de- 
gree of freedom. Each reflection eigenvalue has a twofold 
spin degeneracy. 
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The polar decomposition of the reflection matrix reads 



0.8 




cos OL i sin a\ ( 

. 



I sm a cos a 



(3.14) 



as required by unitarity and the particle-hole symmetry 
condition (j2.5|) . Here U,V are unitary M x M matri- 
ces and ck ~ diag (q;i, a2, . . . aA/) is the diagonal ma- 
trix of real angles that determine the reflection eigenval- 
ues Rn = sin^Q!„. As before, we might have replaced 
sina„ I— >■ \JRn and cosq;„ h- -^/l 



Rn in this polar de- 
composition for class C, but not when we additionally 
impose time-reversal symmetry (in class CI). 

The factor i in Eq. p.l4p accounts for the 7r/2 phase 
shift associated with Andreev reflection of an electron 
into a hole from the opposite spin band. No such factor 
appears in the polar decomposition p.ip in the absence of 
spin-rotation symmetry, because there it can be absorbed 
in the unitary matrices (which in that case contain both 
spin bands for electrons and holes). 

The probability distribution of the Andreev reflection 
eigenvalues in the CQE was calculated previously by 
Khaymovich et ali^ We flnd 

M 

in agreement with Ref. 20 (up to an evident misprint, 
n,^j instead of n,<j)- 



D. Class CI, ensemble T-CQE 
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FIG. 2: Probability distribution of the conductance in the 
CRE, for channel numbers A'' = 1,2,3,4 and topological 
charges Q — —1 (red solid curves) and Q = -1-1 (blue dashed 
curves). In the lower panel the thick vertical lines indicate a 
delta-function distribution. 



Finally, in class CI we have the additional requirement 
(12. 6p of time-reversal symmetry. The polar decomposi- 
tion becomes 

CU \ / cos a isina\ I 0\ 1G^ 
U*J \isma cosa J y WJ' 

The distribution of the i?„'s in the T-CQE (each dou- 
bly degenerate) is again given most compactly in terms 
of the variables f„ = sina„ G (— 1: 1) (with i?„ = ^,^J. 
We find 

M 

^({u)« n i^^-^^-i- (3.17) 

i<j=l 

As in the T-CRE, there is no repulsion between pairs of 
Andreev reflection eigenvalues in the presence of time- 
reversal symmetry. 



IV. DEPENDENCE ON TOPOLOGICAL 
QUANTUM NUMBER OF THE CONDUCTANCE 
DISTRIBUTION 

A. Broken time-reversal symmetry 

From the probability distributions P{{Rn}) in Sees. 
IIII Al and IIII B I we readily calculate the distribution P{G) 
of the conductance (|2.2I) , in both the topologically trivial 
and nontrivial phases. We first consider the case without 
time-reversal symmetry (class D, ensemble CRE). Using 
Eqs. p.7p ~ p.9p we obtain the distributions for the first 
few channel numbers N = 1, 2, 3, 4. The results are plot- 
ted in Fig. [5J and given by: 

• For A^ = 1, the conductance G/Gq — 1 — Q without 
statistical fiuctuationsi^i^ 

• For N ~ 2, the conductance G/Gq = 2 for Q = — 1 
without statistical fiuctuations; if Q = 1, instead 
G/Go = 4g with P{g) = 1. 
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p = l 


p = 2 


(((G 
p = 3 


/Gon 

p — A 


p = 5 


N = 1 


: 2 


: 


: 


: 


: 


TV = 2 


2 : 2 


1-0 


: 


— — 

15 ^ 


: 


N = 3 


3 : 3 


3 • 3 
5 • 5 


2 I 2 
5 ■ 5 


6 • 6 
175 ■ 175 


24 '■ 24 
35 • 35 


N = 4 


4 : 4 


4 : 4 

7 ■ 7 


: 


176 '. 32 
735 • 147 


: 


N = 5 


5 : 5 


5 : 5 
9 • 9 


: 


10 : 10 
2079 ■ 2079 


8 '. 8 

63 ■ 63 



TABLE II: First five cumulants (p < 5) of conductance in the 
CRE, calculated for number of modes N < 5 and topologi- 
cal quantum number Q. (The first entry in each cell is for 
Q — 1, the second entry is for Q — —1.) The conductance 
distribution depends on Q starting from the A''-th cumulant 
(bold). 



• For = 3 the conductance G/Gq = 1 — Q + Ag, 
with P{g) = 3(i - - g f. 

• For = 4 the conductance G/Go = 2 + 4g if Q = 
-1, with P{g) = 3032(1 _ 5)2^ while if Q = 1 we 
have G/Go = 8g with P{g) = 6(1 - |1 - 2g\)^. 

In these expressions, g denotes a stochastic variable in 
the range (0, 1). 

From Fig. [5] we see that upon increasing N, the con- 
ductance distributions for Q — 1 and Q = —1 become 
more and more similar. To quantify the difference, we 
list in Table in the first few cumulants ((Gp)) of P{G) for 
several values of A''. Inspection of the table brings us to 
propose that: 

The cumulant of order p of the N-mode conduc- 
tance in the CRE is independent of the topological 
charge for p < N. 

A proof for arbitrary A'^ is given in App. [BJ 

B. Preserved time-reversal symmetry 

In the presence of time-reversal symmetry (class Dili, 
ensemble T-CRE) we can similarly calculate the conduc- 
tance distribution from Eqs. p. 121) and p.l3|) . For small 
A'^ we find: 

• For N = 2, the conductance G/Gq = 4 without 
statistical fluctuations if Q = —1, while if Q = 1 
we have G/Gq = Ag with P{g) = ig^^/^, 

• For A^ = 4 the conductance G/Gq = 4 -I- 4^ if 
Q = -1, with P{g) = i|(l - 5)^5-'/'; if Q = 1, 
instead G/Gq = 8g with P{g) plotted in Fig. H 
upper panel. (The analytic expression is lengthy.) 

The analogous theorem for the Q-independence of low- 
order cumulants now reads: 




I ' ' ' 1 

2 4 6 8 

G/ Gq 

FIG. 3: Same as Fig. H for the T-CRE. 

The cumulant of order p of the N-mode conduc- 
tance in the T-CRE is independent of the topologi- 
cal charge for p < N/2. 

A proof for arbitrary (even) A^ is given also in App. [BJ 



C. Weak localization and UCF 

Weak localization and universal conductance fluctu- 
ations (UCF) refer to the average and to the variance 
of the conductance in the large-A^ limit. Since the de- 
pendence on the topological charge is nonpertubative in 
A'^, these two effects cannot contain any information on 
whether the superconductor is in a topological phase or 
not. As a check, we have calculated the average (G) and 
the variance VarG = (G^) - (G)^ for TV > 1, directly 
from the probability distribution of the Andreev reflec- 
tion eigenvalues. This calculation also allows us to verify 
a conjecture from Ref. on the UCF in the presence of 
time-reversal symmetry. Since the calculation follows es- 
tablished methods in random-matrix theory)^ we only 
give the results. 

The weak- localization correction SG = G — TVGo to 
the conductance vanishes in the CRE and CQE, while 
(5G/Go = i, -1 in the T-CRE and T-CQE, respectively. 
The UCF are given by VarG/Go = ^,1,2,4 in the 
CRE, T-CRE, CQE, and T-CQE, respectively. These 
Q-independent results are in full agreement with Ref. 0. 

All these results assume that the proximity to the su- 
perconductor does not induce an excitation gap in the 
quantum dot. In the CRE and CQE this is realized by 
the pair-breaking magnetic field. In the T-CRE and T- 
CQE we need a 7r-junction to close the gap: Two NS 
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interfaces, coupled equally well to the quantum dot and 
with a TT phase difference of the superconducting phased 
For a single NS interface in zero magnetic field, the pres- 
ence of an excitation gap does not change the value of 
(5G, but the variance of the conductance is changed into^S 
VarG/Go = 9/4^, with /3 = 1 or ^ = 4 in the presence 
or absence of spin-rotation symmetry. Notice that time- 
reversal symmetry breaking then has only a relatively 
small 10% effect on the UCF^ while in the absence of 
the excitation gap the effect on the variance is a factor 
of two 4 



V. CONCLUSION AND COMPARISON WITH A 
MODEL HAMILTONIAN 

In conclusion, we have shown that the distribution 
P(G) of the electrical conductance in a quantum dot con- 
necting a normal-metal to a superconducting electrode 
has a striking dependence on the topological quantum 
number Q of the superconductor, but only if the number 
of modes N in which the current is injected is sufficiently 
small. In the absence of time-reversal and spin-rotation 
symmetry, the distributions for Q = — 1 and Q = dif- 
fer in the average conductance for iV = 1, in the variance 
for = 2, in the skewness for = 3, and in the kurtosis 
for = 4. More generally, the dependence appears in 
the cumulant of order TV or N/2, depending on whether 
time-reversal symmetry is broken or not. 

The system we have considered (Fig. [1]) is constructed 
to ensure chaotic scattering, which is the requirement for 
application of the circular ensembles of RMT. Systems 
of present experimental focus in the search for Majorana 
bound states have a simpler wire geometry, without the 
quantum dot (Fig.^l inset). Impurity scattering within a 
superconducting coherence length from the NS interface 
can still lead to chaotic dynamics, at least if the number 
of modes is sufficiently small that they are fully mixed 
by the disorder. 

To test the applicability of our RMT results to such 
a system we have performed numerical simulations of 
the model Hamiltonian of Refs. [lOlfTll which describes 
an InAs nanowire on an Al or Nb substrate. The 
Bogoliubov-De Gennes Hamiltonian 



n = 




Hr — Ep 
A* Ep 



A 




Aa 



Ef 



y 

m 



(5.1) 



couples electron and hole excitations near the Fermi en- 
ergy Ep through an s-wave superconducting order pa- 
rameter A. (We have made a unitary transformation to 
ensure that the condition for particle-hole symmetry has 
the form used in the preceding sections.) 

The excitations are confined to a wire of width W in 
the X — y plane of the semiconductor surface inversion 
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FIG. 4: Comparison of the probability distribution of the elec- 
trical conductance as predicted by RMT (dashed curves) and 
as resulting from numerical simulation of the model Hamilto- 
nian (|5.1|) (solid histograms). The simulation is for the dis- 
ordered normal-metal-superconductor junction shown in the 
inset. The number of propagating modes in the normal region 
is A'^ = 2 (lower panel) and N = 3 (upper panel) , while the red 
and blue curves are for topological quantum number Q = ~1 
and Q — +1, respectively. The disorder strength is fixed at 
Uo = 130 Eso for AT = 2 and Uo = 100 E^o for iV = 3. The 
values used for Fermi energy and Zeeman energy (in units of 
Eso) are as follows. For iV = 2: Ep ^ 12, Ez = 3.8 (Q = 1) 
and Ep = 13, Ez ^ 9 {Q ^ -1). For TV = 3: Ep ^ 19, 
Ez = 3.8 (Q = 1) and Ep = 19, Ez ^ S {Q = -1). 



layer, where their dynamics is governed by the Rashba 
Hamiltonian 



P 



U{r) 



2moff ^ ' ' h 

(5.2) 

The spin is coupled to the momentum p = —ihd/dr 
by the Rashba effect, and polarized through the Zee- 
man effect by a magnetic field B parallel to the wire 
(in the x-direction). Characteristic length and energy 
scales are ^so = fi^/iricsaso and Esq — nicsa^^/h^. Typ- 
ical values in InAs are Zgo — 100 nm. Esq = 0.1 meV, 
Ez = ^ffoffMis = 1 meV at B = 1 T. 

We have solved the scattering problem numerically^ 
by discretizing the Hamiltonian (jS.ip on a square lattice 
(lattice constant a — Zso/20), with an electrostatic dis- 
order potential J7(x, y) that varies randomly from site 
to site, distributed uniformly in the interval {—Uo,Uo)- 
The disordered superconducting wire (width VF = 20 a, 
length L — 800 a, A = 4 Ego) is connected at two ends to 
ideal normal-metal leads, obtained by setting A, C/q = 
for X < 0, X > L. The length L was chosen large enough 
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that the transmission probability through the wire was 

< 10-2. 

Results for the probability distribution of the elec- 
trical conductance are shown in Fig. 21 for = 2,3 
and Q = —1,1. (For TV = 1 we simply find the two 
delta-function distributions at G = {e'^/h){\ — Q), as 
expected.) The histograms were obtained by averaging 
over 10^ disorder realizations, conditionally on the value 
of the topological quantum number Q = ±1 (calculated 
from Q = signDet r, as in Ref. [2^) The agreement with 
the predictions from RMT is quite satisfactory. 
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Appendix A: Calculation of the invariant measure 

In this appendix we derive the probability distribu- 
tions of the Andreev reflection eigenvalues in the circular 
ensembles, given in Sec. IIIII We work out the deriva- 
tion for the symmetry classes D and Dili, for an even 
number N = 2M of modes and for topological charge 
(5 = 1, following established methods of random-matrix 
theory;^ The calculations for the other ensembles are 
entirely analogous, so we do not present them here. 

The circular ensembles are characterized by a uni- 
form probability distribution, constrained only by fun- 
damental symmetries. Uniformity in the unitary group 
is defined with respect to the invariant (Haar) measure 
c?/i(r) = r^dr = Sr. Since the polar decompositions in 
Sec. IIIII give a parametrization of the (unitary) reflection 
matrix r in terms of the angles a„ , we can transform the 
measure into d^{r) = J Hi dpi Yin dan- The pi's are the 
degrees of freedom of the matrices of eigenvectors and J 
is the Jacobian of the transformation. From this expres- 
sion the probability distribution of the angles Q!„ follows 
upon integration over the pi's, 

P({a„}) (X J jY[dp,, (Al) 

i 

up to a normalization constant. 

The Jacobian can be found from the metric tensor 17^^, 
which can be read off from the trace Ti Srdr'' when it 
is expressed in terms of the infinitesimals da„ and dpi 
(collectively denoted as da;^): 

Ti-SrSr^ =Y,9t.i^dx^dx,, J = \Betg^,\^/^. (A2) 

We carry out this calculation first for class D and then 
for class DHL 



1. Class D (ensemble CRE) 

In view of the polar decomposition (|3.1|) one has 



\Q V* \ SU* 



lT 'I (A3) 
\ 5V* ^ ' 



where we abbreviated 



The quantities SU — U^dU and 6V — V^dV represent 
measures on the eigenvector manifolds. We used that 
d{V^V) = ^ idV^)V = {SVy = -SV. 
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Substitution of Eq. (jA3l) into Tr Sr6r^ = Tr drdr^ gives 

y SV*J \ SU*J 
- Tr [SU^ + {5U*f + (5V^2 + (5v^*)2] 
+ Trdidit. (A5) 

(All other cross terms vanish.) In terms of F and A this 
can be expressed as a sum of five traces, 

TTiSrSr'f = Tr {T5V - 5UT){T'^ 5U - SVV'^) 

+ Tr {T5V* - 5U*T){T'^5U* ~ 5V*T'^) 
+ Tr {K5V* - SUA){A^dU - SV*A^) 
+ Tr {A6V - SU*A){A^SU* - SVA^) 
+ Tr dLdL'^ 
= Ti+T2 + T3 + Ti + T5. (A6) 

Each of the traces in Eq. (|A6P is of the form Tr AA^ = 

I 

Ti 



We denote A: = 0, 1 for fc 1, 0. 

We group the indices of the matrices SU and SV into 
2x2 blocks, and consider first the off-diagonal blocks. 
For these we can choose as independent parameters 

5U2r,2s, SU2r-l,2s, SU2r,2s~l, <5C^2r-l,2s- 1 > 
SV2r,2s, <5V2r-l,2s, '5V2r,2s-l, <5V2r- l,2s-l , 

with 1 < 7" < s < M. The real and imaginary 
parts, denoted by SU^, SU^, SV^, 5V^, produce a total of 
8M(M — 1) independent parameters. The contribution 



Sij l^u l^i ^-iid is therefore real. Since the second line is 
the complex conjugate of the first line and the fourth line 
is the complex conjugate of the third line, their traces are 
the same, hence Ti — T2 and T3 = T4. For the evalua- 
tion of the expression we need to distinguish between the 
different values of the topological quantum number and 
between the cases of odd and even number of channels. 

We work out the calculation for TV = 2M even and 
Q — 1, when A and F are given by Eqs. p.2p and p.3p . 
The trace T5 is easiest to evaluate, 

M 

T5 = ^ML.,p=45]da?. (A7) 

i,j j=l 

This trace contributes a diagonal block to the metric ten- 
sor and a constant factor to the Jacobian, for M inde- 
pendent real parameters. The other two traces Ti and 
T3 require more work, 



(A8) 



(A9) 



to Tr SrSr'' for given values of r and s has the form 

k=2r-l ;=2s-l 

+ 2a [SV^SU^^-SV^Mi]] 

+ 2b [SV2''r^,,2s-lSU^r,2s + SV2\_,,2s-lSUl^2s] 
+ 2b [SV2l2jU^,_,^2s-l + SVl^2sSUl_,^2s-l] 

where we abbreviated a = — 4 sin sin as and b = 

— 4c0SQ;rC0SQ;s. 



M 1 

= ^ {5(cosar + cosas)2|<5i72.^;,^2s-i-(-l)''+''5T/2r-fe,2«-/l' 

r<s-l k,l=0 

+ i(c0Sa^ - COSasf\6U2r-k,2s-l + h'^)'"^'6V2r^-k.2s-l\^} 
M 1 

+ 51 COS"^ am\SV2rn-k,2rn^l ~ (" 1)''+' ^C^2™-fc,2m-rl ^ 

m k,l=0 
M 1 

= XI {^(sinar + sinas)^|(5V"2*-fc.2s-; - '5C/2r-A;,2s-/P 

r<s = l k,l=Q 

+ ^(sinar - sinas)'^\5V2r-kas-i + ^C^2r-fe,2s-/P} 

M 1 

+ XI Sin^"™l^^2m-fe,2m-i - '5t^2m-fe,2m-/r 

I 
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The contribution to the metric tensor is a block matrix 
with elements 
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-bTQ 
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4to 





Wo 



















4to/ 



where the Pauli matrix and the 2x2 unit matrix tq 
were introduced to account for real and imaginary parts 
in a compact way. The determinant of this matrix is 
(sin^ ar — sin^ oig)^ , hence the contribution to the Jaco- 
bian from the off-diagonal matrix elements is 



Jo 



ff-diasonal 



M 

n 

r<s — 1 



sm Ur 



SIV? Us 



(AlO) 



Next we consider the diagonal blocks. We choose as 
independent parameters 

Wi = -i{5V2m,2m ~ 5U2rn-l,2m-l) , 
W2 — —i{5V2m~l,2rn-l — SU2m,2m), 
W3 = -~iiSV2m,2m + SU2m,2m), 
W4 — <5V2m-l,2m + SU2m,2m-l- 

These are in total 5M real parameters. (Since wi,W2, 
are real numbers they contribute only M parameters 
each.) The contribution to Tr SrSr^ is 

2 2/ 2 \ 2 2 

Wi + W2 + 2['W2W3 — W1W3 — W1W2 + W;^) sin a„i + 2W4, 

and the contribution to the metric tensor is the block 
matrix 



/2 








1 

■ sin^ Q!„ 







sin^ 



. 2 

sm rv 





— sin^ ar, 
sin^ am 



with determinant 2(sinam cosOm)^. Hence the contribu- 
tion to the Jacobian from the diagonal matrix elements 
is 



M 



-'^diagonal = ]^ | siu a„i COS a„ 



(All) 



The number of independent parameters that we have 
accounted for totals to 8Af^ — 2M, which should equal 
the number of degrees of freedom of a matrix in class D. 
The matrix space in class D is isomorphic to the group of 
2A^ X 2N orthogonal matrices;^ which indeed has N{2N~ 
1) — 8AP — 2M degrees of freedom. 

Gathering all terms that contribute to the Jacobian in 
Eq. (|Aip . we obtain the probability distribution 



M M 

sin^ ar - sin^ | | sin am cos | ■ 

r<s— 1 m— 1 

(A12) 

The integration J dpi over the degrees of freedom of 
the eigenvector matrices only contributes a pref actor, 
which can be absorbed in the proportionality constant. 
Upon transformation to the Andreev reflection eigenval- 
ues Rn = sin^ an , we arrive at the result (|3.7p stated in 
the main text. 



2. Class Dili (ensemble T-CRE) 

For the treatment of class Dili it is useful to notice 
the similarity of the polar decomposition of i'SyV given 
in Eq. (|3.10l) to the one of r in class D given in Eq. 
p.ip . Since S{i'Eyr) — Sr all the equations up to Eq. 
(|A9j) derived for class D also hold for class Dili, upon 
replacement U 1-^ ft and F 1— >■ fi*. (As before, we only 
give the detailed derivation for Q = 1.) The expressions 
for the traces Ti and T3 then simplify to 



M 



Ti = 



^ {(cOSttr -I- COSasf [\Sn2r,2s " '5^^2r-l,2s~l P + \Si^2r-l,2s + Sfl2r,2s-l\'^] 
r<s— 1 

+ (cOSttr - COSa^)^ [\5il2r,2s + '5f^2r-l,2s-l P + l'^^^2r— 1,2s - <5f^2r,2s-l P] } 
M 

2^cos2 am\Sn 

2m, 2m m — 1,2m— 1| ; 



M 



Ta = 2 ^ {s'mar - smas)^\Sil2r-k,2 



s-l\ 



r<s=l k,l=0 



r 



(A13) 
(A14) 



For the off-diagonal blocks we choose 

^^^2r,2s, Sil2r-l.,2s, Sil2r,2s-1; '^^^2r- l,2s-l , 



with 1 < 7' < s < M, as the independent real parameters 
(a total of 4M2 - AM). The contribution to Tr SrSr^ for 
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given values of r and s is 

2r 2s 



k=2r~l l=2s~l 

+ 2d [5n. 



^2r-l,2s""2r,2s-l 



""2r-l,2s''"2r,2s-l 



with c = 4 — 4 sin sin and d — A cos cos as . The 
contribution to the metric tensor is a block matrix with 
elements 

/ cro -dT^\ 
ctq dvz 

C?Tz CTo 

V-dr, cro / 

with determinant (sin a,- — sin cks)*. The contribution to 
the Jacobian is 



JoS-d 



iagonal 



M 

n 

r<s— 1 



(A15) 



The diagonal blocks have M independent degrees of 
freedom, 

Si^2m 2m + ^f^2m-l 2m-l, 

which contribute to the Jacobian a factor 

M 



-'^diagonal = W |cOSQ!m| ■ 



(A16) 



The total number of independent parameters (includ- 
ing also the M degrees of freedom from the a„'s) is then 
4M^ — 2M. This agrees with the number of degrees of 
freedom of the matrix space 0{2N)/U{N) in class DIII4 

The distribution of the an's results from the product 

of ^off-diagonal and (/diagonal; 



M 



Isinckr - sinasi'' |cosa,„|. (A17) 



r<s — 1 



Transformation to ^„ — sin a„ gives the expression (I3.12p 
in the main text. 



Appendix B: Proof of the topological-charge 
theorem for the circular ensembles 

The theorem we wish to prove states that the p-ih 
cumulant of the conductance in the A^-mode circular en- 
semble is independent of the topological charge Q for 
p < N/d, with d = 1 in the CRE and d = 2 in the 
T-CRE. 

We start from the definition (|2.2p of the conductance, 
which we rewrite as 



G/Go = iTr[l-rtT,r(l+T,)], 



1 

,0 -1 



(Bl) 



The reflection matrix r is a 2A'" x 2A'" unitary matrix, sat- 
isfying the particle-hole symmetry relation (|2.3p . which 
we rewrite as 



'0 r 



(B2) 



This equation implies that Tr rV^r = 0, hence Eq. (|B1 
reduces to 



G/Go = iTr[l-rtT,rr,)]. 



(B3) 



The p-th cumulant of G contains only averages ruq — 
{[T^T r'^ TzTTzY) with q < p, hence to prove the theorem 
is it sufficient to prove that rrip is independent of Q for 
p < N/d. 

We first do this for the CRE. Then the average can 
be written as 

mp = J dfi{r){TTr^TzrTzY^il + QDetr), (B4) 

where d^i{r) is the invariant measure of class D. The 
defining property of this measure is that dp,{Ur) = 
d^{rU) = d/i(f) for any 2A^ x 2A^ unitary matrix U that 
satisfies U — TxU*Tx- What we seek to prove, therefore, 
is that 

j dfi{r) (TrrV^rr^j^'Detr = if p < A^. (B5) 
We decompose 

N 

r. =E^^"^ ^if = SkASk,n - Sk,n+N) (B6) 
n=l 

and apply this decomposition to one of the r^'s in Eq. 
(lB5l). 



N N N 
Pi=0p2=0 pjv=0 



N 

X 5p^p-^+p2-\ hPN 



Pn 



(B7) 



n=l 

Consider one of the terms 

N 



M= dfi{r) f[ (^TrrV("Vrz)^" Det r. (B8) 

n=l 

li p < N, there is at least one index no £ {1,2,... A^} 
such that pn„ = 0. Transform r i-> [/(""V, with 



6k,i if fc 7^ no,no + A^, 
Sl,io+n if k = no, 
Si,no if fc = no -I- A^, 



(B9) 



a real, symmetric, unitary matrix which commutes with 
Tx- This transformation does not change the invariant 
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measure, dii{U^'^°^r) = dfi{r), while the integrand trans- 
forms to 



M= d^l{r) W (^Trr^C/("'')r(")i7("°) 



X Detf7("''V 

N 



Mr) n (TrrV(" 



rTz 1 Det r 



= -M, 



(BIO) 



since Det C/^""' = —1 and [/("") commutes with r*^"^ for 
n ^ no, while png = 0. Hence = 0. 

This completes the proof for the CRE. For the T-CRE, 
we seek to prove that 

J d^(r) (TrrV^rT^)^PfiEj,r = if p < N/2, (Bll) 

where now dfj,{r) is the invariant measure of class DHL 
The invariance property reads dfi(T,yU'^TjyrU) = dfj,{r) 
for any 2N x 2N unitary matrix U that satisfies U — 
TxU*Tx- Since Tz and Sy commute, we may rewrite Eq. 
pTTj) as 

j d^i{r) {Tvr'^T.yT^T.yrT^Y FiiT^yT = Q if p < N/2. 

(B12) 

Substitute the decomposition (|B6p in both the Tz's, 



N N 



p\ 



Pll=0 Pjvjv 



nn,mP" 



N 



P-J2rz m Pi 



n.m—l 



(B13) 



Consider one of the terms 



r ^ 

M= dfi{r) Yl (Trr1'Sj,r(")Eyrr(™) 



X Pf iSj^r. 



(B14) 



If p < N/2, there is at least one index uq E {1,2,... N} 
such that Pngm — and Pntig = for each n,m G 
{1,2,... TV}. Transform r j:yU'-"'>^T,yrU^"''>\ with 
jjino) defined in Eq. (jB9|) . This transformation does not 
change the invariant measure, so the integral transforms 
to 



Jdn{r) f[ (Trr^EyC/("'''r(")C/("'''Sj,r 

n,m—l 

X [/(no)^(m)^(no)^^"" pf (f7("o)j5]^^^(no)^ 

r ^ 

= - dfiir) Y[ (TrrtSj^r^^^S^rr^™)' 



n,m— 1 



X PfiEj^r = -X, (B15) 

where we have used that PfXYX'^ = (DetX)(Pfy). 
Hence A4 = and we have completed the proof for the 
T-CRE. 
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